Quantum dot spin qubits in Silicon: Multivalley physics 
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Research on Si quantum dot spin qubits is motivated by the long spin coherence times measured 
in Si, yet the orbital spectrum of Si dots is increased as a result of the valley degree of freedom. The 
valley degeneracy may be lifted by the interface potential, which gives rise to a valley-orbit coupling, 
but the latter depends on the detailed structure of the interface and is generally unknown a priori. 
These facts motivate us to provide an extensive study of spin qubits in Si double quantum dots, 
accounting fully for the valley degree of freedom and assuming no prior knowledge of the valley- 
orbit coupling. For single-spin qubits we analyze the spectrum of a multivalley double quantum 
dot, discuss the initialization of one qubit, identify the conditions for the lowest energy two-electron 
states to be a singlet and a triplet well separated from other states, and determine analytical 
expressions for the exchange splitting. For singlet-triplet qubits we analyze the single-dot spectrum 
and initialization process, the double-dot spectrum, singlet-triplet mixing in an inhomogeneous 
magnetic field and the peculiarities of spin blockade in multivalley qubits. We review briefly the 
hyperfine interaction in Si and discuss its role in spin blockade in natural Si, including intravalley 
and intervalley effects. We study the evolution of the double-dot spectrum as a function of magnetic 
field. We address briefly the situation in which the valley-orbit coupling is different in each dot due 
to interface roughness. We propose a new experiment for measuring the valley splitting in a single 
quantum dot. We discuss the possibility of devising other types of qubits in Si QDs, and the size of 
the intervaley coupling due to the Coulomb interaction. 



I. INTRODUCTION 

A practical QC architecture requires scalability and 
long coherence times. Spin-based solid state quantum 
bits (qubits) are known to have long coherence times,^^ 
while also offering the promise of scalability, and are nat- 
ural building blocks for quantum computation. The con- 
ventional criterion establishing an acceptable coherence 
time for fault-tolerant quantum computing is for it to al- 
low lO"* operations. With the operation time set by the 
gating time « Ins the coherence time required is « 10/iS. 

Phosphorus donor nuclei in Silicon have been known 
since the 1950s to have some of the best spin coherence 
properties in solidSj'SHS! and were therefore proposed as 
viable candidates for qubits.^ However, fabrication of or- 
dered and gated donor arrays and coherent control over 
donor electrons has turned out to be extremely difhcult 
(note nevertheless outstanding recent developments.^"^") 

Electron spin qubits in quantum dots (QDs), ^partic- 
ularly in GaAs, have been actively studied in the past 
decade, and c ohere nce times in excess of « 100/xs have 
been achieved .'^'^lill Various schemes exist for implement- 
ing spin qubits in QDs. For the original single-spin 
qubits, one electron spin is initialized in each dot in a 
QD arrayjf^ with spin up and down representing the two 
states of a qubit. One-qubit rotations are accomplished 
using magnetic resonance techniques and two-qubit rota- 
tions by means of a pulsed exchange interaction. One can 
also use selected two-spin states to encode a logical qubit: 



in singlet-triplet qubits it is the singlet and unpolarized 
triplet states that are used for encoding. Here exchange 
coupling gives t he spli tting of the two-spin states mak- 
ing up the qubit .'^'^Sllll One-qubit rotations are currently 
carried out using an inhomogeneous magnetic field across 
the two dots and pulsed exchange splitting. Two-qubit 
operations use the different dipole moments of the singlet 
and triplet states. Larger encoding schemes such as us- 
ing three-spin states have also been proposedj^i^ aiming 
at an all-electrically-controlled architecture. 

Recent years have seen significant experimental 
progress involving single-spin properties s uch a s co- 
herent control, coherence, and measurement pMUHH 
GaAs double quantum dots (DQDs) spin blockadtP 
and charge sensors (quantum point contacts and radio- 
frequency quantum point co ntacts^P^ enable observation 
of single/two-spin dynamics .I^^^IIIISIIII Spin coherence in 
GaAs QDs is mostly limited by hyperfine interaction with 
the nuclear spins,'2SH2Sl though proposals exist for manip- 
ulating nuclear spin states to extend ele ctron spin coher- 
ence timeSj'^^SlI and recent experiments^'^-' have shown 
that dynamical decoupling schemes^^ 39 - 42 (,q^\^ help 
unwind the nuclear spin dynamics. Progress has also 
been made on controlling spins in GaAs QDs by optical 
means, such as the preservation of spin coherence, ^^^^ con- 
trol of the nuclear spin polarization,S3Hifil and ultrafast 
spin rotations.'^ Yet experiment is far from achieving re- 
liable control over the nuclear field in GaAs. 

Silicon has outstanding spin coherence propertieJ^] due 
to small spin-orbit coupling,'** small hyperfine interaction 
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with nuclear spin^^ (which can be reduced by isotopic 
purification^) and absence of piezoelectric coupling!^ 
The spin coherence time T2 (after Hahn spin echo) for 
donor electron spins in bulk Si:P has been reported to 
range from 60 ms^ to 300 ms.^°'^ This is the longest coher- 
ence time measured in electron spin qubits, and greatly 
exceeds the values reported in GaAs QDs, which, after 
Hahn echo, range from ~ to ~ SO/zs.-^ Isotopic 

purification of Si, where ^^Si nuclei are preferentially re- 
moved in favor of ^^Si and ^"^Si nuclei from natural Si . al- 
lows enhancement of Si electron spin coherence timesP^ 
The maturity and continuous innovation in Si microfab- 
rication could be of great help in scaling up a Si-based 
QC architecture. Work is under way on Si QC architec- 
tures b ased on donor electron or nuclear spins in a Si:P 
system,'^'^ single electron spins in gate-defined quan- 
tum dots in Si/SiGtP^mi or Si/SiOsP^ and holes in 
SiGe nanowires.'^SI Experimental progress has been made 
in b oth donor-bas ed devices^ ^*' and gate-defined 

QP3 l8llOI21l56l57l66l t68l 

The valley degree of freedom presents a potential ob- 
stacle to spin QC in Si. In bulk Si there are six de- 
generate conduction band minima near the X points in 
the first Brillouin zone. Impurity scattering breaks the 
symmetry of the lattice and allows coupling between the 
valleys, and in samples in which two or more donor cen- 
ters are present interference becomes possible and is sen- 
sitively dependent on the exact atomic locations of the 
two donors.'S^^ Confinement in the ^-direction and/or 
uniaxial tensile strain at the interface cause the valleys 
perpendicular to the interface to have a smaller energy 
by several tens of meV than the valleys in the plane of 
the interface. A sharp interface further couples the two z- 
valleys by producing a valley-orbit coupling A (a complex 
number) .^2! A is generally not known a 'priori^ is sample- 
dependent,'^ and is different for different interfaces, such 
as Si/SiGe an d Si/S i02. Various measurements of A have 
been reported^^^^ and efforts have been devoted to 
calculating A.ElE2lI3IHi] Yet the fact that A is not known 
and no standard experiment exists to measure it in a QD 
at low field is a strong motivating factor for our work. 

In a recent articl^^l .^^g discussed singlet-triplet qubits 
in Si QDs, emphasizing the role of the valley degree of 
freedom. The focus of Ref.l82lwas on determining the fea- 
sibility of an experiment analogous to Ref. |T]in a Si/Si02 
(or Si/SiGe?^ DQD. It was found that the valley degree 
of freedom affects qubit initialization, operation and spin 
blockade. For |A| >• ksT a singlet-triplet qubit can be 
implemented in the same way as in GaAs. For | A| « fc^T 
initialization will be problematic, since, unlike in GaAs, 
the energy scale is no longer set by the confinement en- 
ergy but by the valley splitting (we note a similar dis- 
cuss ion o f spin blockade in carbon nanotube quantum 
dots.™l) 

In this work we present a detailed analysis of both 
single-spin and singlet-triplet qubits in multivalley Si 
QDs within the effective mass approximation. The dis- 
cussion of valley physics is general and applies to all 
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Figure 1: (color online) Conduction-band valleys in silicon 
in the bulk and in 2D. (a) In bulk Si there are six equivalent 
conduction-band minima positioned along the Cartesian axes, 
(b) Uniaxial strain at the interface causes the valleys paral- 
lel to the interface (i.e. in the xj/-plane here) to be higher 
by several tens of meV than the valleys perpendicular to the 
interface. The interface potential further splits the lower- 
energy valleys by an amount A, which present experimental 
work is striving to determine unambiguously. Based on mea- 
surements to date we expect A <lmeV. 



types of Si QDs, yet when concrete figures are needed 
we focus on Si/Si02. Many aspects of the following con- 
siderations apply also to carbon-based QDs, where the 
valley degree of freedom is present as welLl^^HlIl "We con- 
sider first single-spin qubits, i.e. the Loss-DiVincenzo 
(LDV) architecture.'^ We determine the spectrum and 
the conditions required for implementing single-spin 
qubits and coupling adjacent qubits by means of the ex- 
change iriteraction. We subsequently study singlet-triplet 
qubits.UJi^ Part of the discussion is intended as a fol- 
low up to Ref. 'S?, providing further theoretical details 
of singlet-triplet qubits including all the relevant ana- 
lytical formulas and matrix elements. The remainder 
contains an extensive discussion of the role of the hy- 
perfine interaction in spin qubits in natural Si, includ- 
ing intervalley effects, an exhaustive discussion of single- 
qubit gates and scalability, as well as new experimental 
initiatives. We conclude that in both schemes, for any 
A, sm(/le-qubit operations are feasible. If |A| < k^T in 
single-spin qubits exchange-based entanglement will not 



3 



work since exchange between pairs of states with differ- 
ent vaUey composition is vastly different. For this reason, 
the proposed Coulomb interaction-based two-qubit oper- 
ations in the singlet-triplet scheme will not work either. 
The outline of this article is as follows. In Sec. |TT] we 



introduce the DQD Hamiltonian. In Sec. Ill we study the 
effect of the valley degree of freedom on single-spin qubits 
and in Sec. |IV| we investigate singlet-triplet qubits and 
the role of the valley degree of freedom in initialization, 
coherent manipulations, and readout. We give a brief 
review of hyperfine interactions in Si, then discuss in- 
travalley and intervalley effects relevant to singlet-triplet 
qubits. Interface roughness is discussed in Sec. |Vj other 
types of qubits in Sec. VI and in Sec. VII an experi- 



ment is discussed that can estimate the valley splitting. 
We also discuss briefly the effect on spin qubits of inter- 
face roughness, which may lead to a different valley-orbit 
coupling in the two dots. We end with a summary and 
conclusions. 



II. MODEL OF A DOUBLE QUANTUM DOT 

We focus on a DQD in a two-dimensional electron 
gas (2DEG) in the inversion layer of a MOSFET grown 
along the z-direction (which is also the crystalline (001) 
direction). The dots are located at Rd = {X£),0,0), 
where D — L,R, and Xn = Xq while Xl — —Xq 
(since the z and —z valleys are isotropic in the xy 
plane, we can choose the x and y axes by convenience). 
We do not consider spin-orbit interactions. The DQD 

Hamiltonian H = Hq + J2i=i2^^^^^ where rep- 
resents the valley-orbit coupling discussed below and 
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Vee, where T is the kinetic 



^0 = E.= l,2 {T^' 

energy operator and Vq the confinement potential 
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-eEx, 
(1) 

with a and h the Fock-Darwin radii for in- and out- 
of-plane confinement respectively. The external electric 
field E II i; is used to bias the DQD and is a crucial tuning 
knob for singlet-triplet qubits. The Coulomb potential 
between electrons at ri and r2 is Vee = e^/(eki — '"2I), 
where the effective dielectric constant e = (esi + esj:o2)/2 
accounts for the extra screening due to the image charge 
in Si02 (for Si/SiCe such averaging is not needed). We 
use a quadratic potential to model the growth-direction 
confinement, which is generally triangular. The in-plane 
confinement potential for each dot is harmonic near the 
center, Vd{x) — {hujo/2a'^) {x ~ Xu)'^ , with Gaussian 
ground states. 

The confinement along the growth direction splits the 
six-fold degenerate Si conduction band into a doubly- 
degenerate branch of lower energy and a fourfold- 
degenerate branch of higher energy. The two lower en- 
ergy valleys ^ = z, z are along the direction perpendicular 
to the interface at fe^ = ztko z, where kg = 0.85{2tt /asi), 



and the lattice constant asi = 5.43A. The ground-state 
wave functions satisfy (T + Vd) — Eq D^, with 

D^^FDir-RDje'^^-^u^ir), (2) 

where D ~ L 01 R. The envelope functions are 
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7r3/4(a2&)i/: 



(3) 



where a = 



is the in-plane Fock-Darwin radius, 5, 



the growth-direction confinement length, = 0.191too 
is the in-plane (transverse) Si effective mass, and mg the 
bare electron mass. The lattice-periodic function (r) = 
c\^&^^'^ ^ with K the reciprocal lattice vectors. 
Overlaps of the form {L^\R^) are given by / = e~^^, 
where d = Xo/a. Overlaps of the form {D^\D^^) are 

suppressed by an exponential of the form e ^ ^ 1, 
where Qz = — 2fco (n^ is an integer) . With h w 3nm 
typically,'2£! h'^Ql/A « 800 when = and « 150 
when Uz — \. Matrix elements of potentials smooth in 
space and interactions between states from different val- 
leys are suppressed by such prefactors: matrix elements 

of the form {D^\Hq\D_^) and (D^^^D^ ^ |7?o|£'i^££'f 
The one exception is the valley-exchange Coulonib inte- 
gral below. Two-particle matrix elements in which 
each electron belongs to a different valley, for example 

(Z?l^''iI'^^|i?ol^^^''^i^^)j are of the same order of magni- 
tude as tnose in which all functions belong to the same 



valley, such as {D^^ Df\HQ\D'l'' D'l^') . Matrix elements 
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of interactions on the atomic length scale are not ex- 
ponentially suppressed. One prominent example is the 
hyperfine (hf) interaction between a conduction electron 
and a ^^Si nuclear spin (in a natural Silicon substrate 
about 4.7% of the nuclei are ^^Si). 

Based on multiple experiments conducted over the past 
50 years and on the recent calculations (see Sec. IV B 4 



for details), it is assumed that the hf interaction in Si is 
dominantly of the contact character, i.e. the Hamiltonian 
is given by 



Hhf = AqVq ^ S • I,<5(r - R,) , 



(4) 



where the sum runs over the positions of the ^^Si nuclei, 
= |A'ofi.^7s7j/i'o (with 75 and 7,7 being the electron 
and nuclear gyromagnetic factor, respectively), vq is the 
unit cell volume and 1^ are the nuclear spins. The valley- 
diagonal matrix element is 



{D^\HHFm 



(5) 



where the hf interaction energy is = AqT] with 77 = 
|u^(Ri)P ~ 160, see Ref. |49] and references therein 
(^ ~ 2 /ieV in Si) . This Hamiltonian provides a valley- 
orbit coupling mechanism, i.e. we have {D(^\Hhf\D-() = 
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AoiyQ J2i (R-i)-C'-^(Ri)S • li- The consequences of this 
couphng are explored below. 

The interface potential couples the z and z valleys if it 
is sufficiently sharp on the atomic scale. This is contained 
in Hy , which represe nts the valley-orbit coupling A stud- 
ied in detail recently.'^^'^ The magnitude and phase of 
A are of great interest in Si research. It has been shown 
that | A| can range from zero to 0.25 meV (maybe even 
larger ),'221Sil and both its magnitude and phase can be 
controlled by an electric field. In the present work 
is a single-particle phenomenological coupling between 
the valleys {D^\H^\D_^) = A = |A| e"''-^, with |A| > 0. 
Throughout this work hwz ^ hujo ^ |A|. 

We diagonalize the single-particle Hamiltonians in- 
cluding Hy and obtain the valley eigenstates D± — 
(l/\/2) {Dz ± e"l'Dg) with eigenvalues Eq ± |A| in each 
of the quantum dots. We construct orthogonal (Wan- 
nier) single-dot wave-function^SH j^^ = Li-gRti ^ 

Pl ^ Rj-gLj = (1 - VT^)/l, so that 

{Ri\L^) = 0. We define 

£0= {D^\T + Vd\D^) 

A= {D^\Hy\D^^) (6) 
{L^\Ho\R^) + {L^L^\V,,\L^R^). 

Orthogonalizing L± and R± we obtain L± — 9R± ^ 
and R+ — -^±^g-^± ^_ These are the single-electron 

y/l-2lg+g^ & 

states we will use henceforth (e.g. to construct two- 
electron states.) The orbital excitation energy is assumed 
large (> ImeV, or > |A|.) 

III. SINGLE-SPIN QUBITS 

For single-spin qubits we first consider initialization. In 
single- valley systems the lowest energy level in zero mag- 
netic field is spin-degenerate but orbital-non-degenerate. 
Initialization of a definite single-spin state requires a 
magnetic field to split the spin degeneracy via the Zee- 
man interaction, so that selective tunnelling from a reser- 
voir can have high fidelity. Experimentally in GaAs this 
process requires magnetic fields of several Tesla due to 
the inherent smallness of the Bohr magneton, the small 
g-factor in GaAs {\g\ ^ 0.4), and the constraint of ^ 100 
mK electron temperature in the reservoir. In Si g ^ 2, so 
that at dilution refrigerator temperatures much smaller 
magnetic fields will be required than in GaAs, rendering 
spin initialization a relatively simple task. 

In Si however there are two possibly closely-spaced one- 
electron levels D+ and D_, separated by 2|A|. Initializ- 
ing into a definite orbital state is dependent on the size of 
I A|. Since the interaction of spins in states and D_ 
with a magnetic field is described by identical Zeeman 



Hamiltonians, single-qubit rotations can be performed no 
matter which state is initialized. 



A. Two-qubit gates 

In the LDV architecture two-qubit gates are imple- 
mented using the exchange interaction between neigh- 
boring dots. Two-electron wave functions are superpo- 
sitions of spin singlet and triplet states. Singlet wave 
function j^g) — |S')|xs), where 15) is a two-particle spa- 
tially symmetric function and \xs) is the two-spin sin- 
glet state given by |xs) = (1/v^) | t^^U^^^ - I^^H^^^ 
). The superscript (i) here denotes the i-th electron. 
Triplet wave functions |^t) — \T)\xt), where |T) is a 
two-particle spatially antisymmetric function, and \xt) 
represents the spin-triplet states, and is an abbrevia- 
tion for |x+) = I t^'H^'^), IX-) = I l^'U^'^), and 
Ixo) = (1/^2)1 t^'U*'^ + ;('^t^'^)- The wave func- 
tions |x±) represent the polarized spin triplets and are 
given by | f (i)f (2)^ a,nd | i'-^^i^'^^) respectively, while |xo) 
represents the unpolarized spin triplet and has the form 
Ixo) = (l/%/2)| f^h^^^ + ;(i)t^^^)- We will frequently 
refer to just \S) and |T) as the singlet and triplet states 
respectively, because the orbital part is where valleys play 
an important role. 

Translating the LDV architecture to multi-valley sys- 
tems requires taking into account the valley degree of 
freedom in the construction of many-particle wave func- 
tions. In general, for a given number N of orbital states 
one can form N{N + l)/2 singlet and 3N{N -l)/2 triplet 
two-electron states. In a single-valley system the num- 
ber of orbital states is = 1, resulting in the trivial 
case of one singlet. In a single Si QD we have N — 2. 
For two electrons we have 2(2-|-l)/2 — 3 singlet states 
and 3[2(2-l)]/2] = 3 triplet states. In a Si double dot, 
when considering the lowest orbital state in each valley, 
N = 4, so that there are 4(4-|-l)/2 = 10 singlet states 
and 3[4(4-l)]/2 — 18 triplet states. Therefore there are 
altogether 28 states for the double dot, of which 10 are 
spatially symmetric (spin singlets) and 18 are spatially 
antisymmetric (spin triplets). The total number of 28 
states corresponds to the combination Cf, where 8 = 2^ 
arises as the product of a factor of 2 for the spin degree of 
freedom, a factor of 2 for the valley degree of freedom and 
a factor of 2 for the two dots. The combination Cf is thus 
the number of combinations to construct a two-electron 
wave function out of 8 possible states, and represents 
the dimensionality of the representation of the permuta- 
tion group on the combined two-dot two-spin two-valley 
SU(2)xSU(2)xSU(2) Hilbert space. 

The singlet and triplet states in a Si DQD can be 
divided into three branches: states with the same val- 
ley composition, comprising wave functions of the form 
{++) and ( ), and states with mixed valley composi- 
tion, comprising wave functions of the form (-1 — ). There 
are thus in total six uncoupled branches, which can 
be labeled {++)S, {+-)S, ( )S and (++)f , i+-)f, 
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Table I: Number of singlet and triplet states in difltoront 
branches of the spectrum in LDV and ST schemes. 





++ 


+- 




Total 


LDV singlet states 


3 


4 


3 


10 


LDV triplet states 


3 


12 


3 


18 


ST singlet states 


2 


3 


2 


7 


ST triplet states 


3 


9 


3 


15 



( — )f (see Table The spatial parts of the ±± singlet 
and triplet states are 



qLR _ 

^±± — 



1 

71, 



(7) 



1 



(1) £.(2) _ £(2)^(1)^ 



For the H — states, the spatial parts of the singlets are 



5^«^_L (2(1)^(2)^^(2)^(1)^ 



1 



(1)r(2) 



+ /2 V 



1 

71 



■^(1)^(2) 



(8) 



form 



H±± = 2eo ± 2|A| 



fk+j iV2 iV2 \ 
iV2 u j 
iV2 j ii 

V k-jj 



All terms are given in the Appendix. In the 



h1'^ = 2eo 



/ k ±j i i \ 
±j k ±i ±i 
i ±i u ±j 
V i ±i ±j u j 



(10) 
branch, 



(11) 



where + applies to and — to H^_, and is 

in the basis {S^'^,S^l,S^?:,Slt} while for Hl_ re- 
place S T. Since the overlap between states from 
different valleys is negligible, matrix elements of the form 

{S+-\Hi^\S^ ) and {T^ I^^oI^h ) are equal up to a sign 

(we will discuss below the case when this condition is not 
satisfied.) The eigenvalues of and are exactly 
the same. In other words, there is no exchange splitting 
between the singlet and triplet states if the two electrons 
are in different valley eigenstates. 



B. Scalability 



while the spatial parts of the triplet states are 
1 

71 



T: 



LR 



1 



)(i)n(2) 



i(2)n(i)^ 



(9) 



V2 



Because we are using the Wannier states D±, there are 
no overlaps such as {L ± \R±). The spatially symmet- 
ric functions S split into three single-valley states from 

the lower-energy branch, three single-valley states 

from the higher-energy H — h branch, and four states from 
the H — branch. The spatially antisymmetric functions 
T split into one single-valley state from the lower-energy 
branch, another single-valley state from the higher- 
energy -|--|- branch and four states from the H — branch. 
As we discussed before, the smooth quantum dot confine- 
ment potential and the electron-electron Coulomb inter- 
action do not introduce further transitions between val- 
leys (the matrix elements are mostly exponentially sup- 
pressed, and the only Coulomb matrix element that is 
not exponentially suppressed is nonetheless very small, 

as we show in the Appendix). Thus the , and 

H — branches do not mix, allowing us to analyze them 
separately. 

The two-electron Hamiltonians in the + + / — — 
branches, in the bases {S^^, S^^, S^^,f^^}, have the 



The structure of the two-electron Hamiltonian has im- 
portant consequences for the operation of spin qubits 

in silicon. The -I— I- and branches have the same 

structure as the single-valley Hamiltonian except they 
are shifted up and down by 2|A| respectively. These 
branches can be regarded as replic as of the usual single- 
valley Hund-Mulliken Hamiltonian.l^SIlZ! The singlet and 
triplet states have the same form as in the single-valley 
case. The principal difference from the one-valley situa- 
tion emerges when we consider the H — branches, where 
we have four sets of degenerate singlet and triplet states. 
This degeneracy is the direct manifestation of the fact 
that exchange integrals for electrons from different val- 
leys vanish (either exponentially small or very small due 
to Bloch function symmetry, see Appendix). 

The essence of exchange gates in the Loss-DiVincenzo 
architecture is the control of the phase difference (intro- 
duced by a finite exchange splitting) between singlet and 
triplet states. In the or -I— I- branches the two elec- 
trons are in the same valley eigenstate, so that we exactly 
recover the single-valley physics of GaAs, and exchange 
gates can be implemented. However, in the -I — and — h 
branches singlet and triplet states are degenerate, ren- 
dering exchange gates impossible. Now assume that the 
two electrons were initialized into two initially indepen- 
dent quantum dots with an arbitrary combination of val- 
ley eigenstates. When the inter-dot barrier is lowered, 
if |A| is sufficiently small, the two-electron states would 
generally populate all four branches (++, H — , — h, and 
). Since exchange gates cannot be implemented in 
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the H — and — h branches, they cannot be implemented 
for these two electrons in general. The only way to per- 
form exchange gates is to have the two-electron states 

populating only the -|--|- and branches. With the 

electrons in a reservoir in a structure-less Fermi distri- 
bution, the only single-electron state that can be loaded 
exclusively is the — valley eigenstate (if 2|A| ^ kgT)^ 
so that when inter-dot tunnelling is turned on, the two 
electrons are in the branch. This is the only vi- 
able option where exchange gates can still be performed. 
The exchange splitting between the lowest-lying singlet 
and triplet states in the — — branch is given by the 
Hund-Mulliken expression j ^^l^^l which can be obtained by 
straightforward diagonalization of Eq. (10), 
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J 



(u-k) 



+ i \l{u-kY + l&P. (12) 



In the H — h and branches there are four energy levels: 

a singlet, a triplet which is higher by J, and two sin- 
glets which are much higher in energy due to the on-site 
Coulomb interaction. In the H — branch one has the same 
four levels (bar the valley splitting), except that for each 
energy there is a singlet and a triplet. In this branch ex- 
change is zero, and the next highest state is separated by 
J. Therefore exchange in the three branches is either J 
or zero. This fact becomes important when |A| < kgT. 
Firstly, for A = 0, the lowest state is a (sixfold) de- 
generate singlet /triplet. For A nonzero but |A| < fc^T, 
one may in general initialize a superposition of states 
and exchange is not well defined. Consider two qubits 
{aLL+ + hL-l and (afli?+ + (3rR.). If ul = = 1 
the states are L+ and R-, and the exchange coupling 
between them is zero. For arbitrary ao, (3d exchange 
depends on the specific values of and /3d and cannot 
be controlled, precluding QC operations. 

Our analysis shows that in the LDV architecture for 
multivalley systems single-qubit operations are straight- 
forward, while exchange-based two-qubit operations are 
not feasible unless |A| ^ fcsT. The optimal strategy 
is to strive to obtain a large valley splitting and work 
towards replicating the single- valley situation. 



IV. SINGLET-TRIPLET QUBITS 

We have studied the implications of multiple valleys 
on the implementation of singlet-triplet qubits in Ref. [521 
In this Section, aside from reviewing initialization (Sec. 
IV A ) , we provide the important technical details and ex- 



plore new features such as the evolution of the spectrum 
as a function of magnetic field, intravalley and interval- 
ley effects due to the hyperfine interaction, the role of 
these terms in the dynamics of singlet-triplet qubits, the 
difference between applied and hyperfine magnetic fields, 
spin blockade and two-qubit operations. 

We use the notation (n^m) to indicate the occupancy of 
the left and right dots respectively. The single-particle 



Figure 2: Two-electron energy levels on a single dot at fi- 
nite A and magnetic field such that 2|A| > Ez as introduced 
in Ref. 1821 In this case the lowest energy state is the sin- 
glet S , followed by the triplet , the degenerate singlet 

5+ -/triplet T]^! and triplet Tjii, and finally the singlet S++. 
Electron spin orientations in T are indicated by arrows. 



energies on the left and right dots are different, — 
{R^\{T+Vq)\R^) and El = {L^\{T+Vq)\L^} . The dimen- 
sionless detuning is defined as 5 = {el — ER)/{2dEo), and 
the dimensionless critical detuning Sc = {u — k)/{2d£Q). 



A. Initialization 

Our study of the singlet-triplet qubit follows the suc- 
cessful experimental implementation in GaAs DQDs, 
which has 4 basic ingredients: initialization in a single 
dot in the (0,2) regime, single-qubit manipulation in the 
(1,1) regime, two-qubit manipulation via electrostatic in- 
teraction, and spin measurement via spin blockade.^ 

First a (0,2) state is initialized {E > 0.) Using the 
notation of Sec. |III[ the four lowest-energy two-particle 
spatial wave functions are. 



oRR _ n(2) 
— -'''ill -^ih 

= (l/N/2)(i?«i?L'^ + i?f ^-'') 



(13) 



ff^ = (1/^2) (i?V^i? 



(1)6(2) 



In the basis {S^^, S^^, Tf/^, 5ff } the Hamihonian is 



Hi 



dot 



2en 



/-2|A| \ 



V 2|Aiy 



(14) 



For A = and no external magnetic field B, all six 
lowest-energy levels are degenerate and it is not possible 
to load any particular two-electron state. An external 
magnetic field splits the triplet energy levels by the Zee- 
man energy Ez (Fig. [2]) Initialization of a (0,2) state re- 
quires an outside reservoir with a Fermi energy Ep (ther- 
mally broadened by sa ksT) tuned to be on resonance 
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with S^^. If |A| > ksT (at T=100mK we require A « 
O.lmeV) then S!^^ can be loaded exclusively. 

The qubit of Refs. |1|4| is in a single-valley system, so 
R± — R. Only a singlet can be made out of the lowest 
lying wave function. The triplet involves higher orbitals 
and is separated by Eq » ksT, minus a small exchange 
term. In GaAs this energy scale is « 1 meV, which is 
much greater than the experimental thermal energy of 
0.01 meV (lOOmK). In Si the energy scale is set by A. If 
A is sufficiently large the physics is that of a single- valley 
problem. 



In the evaluation of Eq. ( 14 1 we have also encountered 



a valley-exchange Coulomb integral 



(15) 



This integral is not suppressed by an exponentially small 
prefactor, yet it is ^ l/icV and we do not take it into 
account. Its evaluation is discussed in Appendix [P} 



B. One-qubit manipulation 

Single-qubit operations on the singlet-triplet qubit are 
performed in the (1,1) regime.'^ The two-particle states 
are the same as those enumerated in the case of single- 
spin qubits, except here we do not include the high- 
energy (2,0) states of the form LL (see Table |l]) In 
the bases {S^^,S^'^,f^^} the matrix elements of the 
Hamiltonian in the + + / branches are given by 

/ j 

ff±±=£L+£i? + ^±2|A|+ ^^/2 ~2dea{5 - 5^) 

V ~h 

(16) 

The matrix elements in the H — branch, in the bases 
^f«} and {ff«,f^f,ffi^}, are 



S,T 



±j t 
±j ±i 

1 ±i -2deo(<5- 



(17) 



where the ± signs apply to the singlet/triplet Hamilto- 
nians respectively. 

The H — singlet and triplet branches always yield the 
same energies. Moreover, for A = the lowest en- 
ergy state is sixfold degenerate: the lowest energy ++, 

H — and singlets, plus the lowest energy H — triplet. 

When the two electrons are in the same valley, the singlet 
and triplet energies differ by J. 

For clarity we focus on a concrete example, taking a Si 
DQD with a=8.2nm, b=3nm, d=2.45 and |A|=0.1meV. 
The energy levels of the system with these parameters are 
plotted in Fig. [s] (analogous to Ref. |82| as a function of 
the dimensionless detuning S. At low detuning there are 
four (0,2) high-energy levels, indicated by the two solid 
lines (representing singlets of the form _R+_R+, _R_i?_) 



and one dashed line (representing two degenerate singlet 
and triplet of the form i?+i?_). The separation of these 
levels is 2|A|. The lower-energy (1,1) levels are a degen- 
erate singlet/triplet of the form I/+-R+ (top solid line), 
a degenerate singlet/triplet of the form L_i?_ (bottom 
solid line) , and two degenerate valley mixing singlets and 
triplets of the form L_|_i?_ and L^R^. For the param- 
eters considered in this example i « 0.02meV, giving a 
splitting at the avoided crossing of « 0.06meV. 



1. Double dot spectra in different parameter regimes 

Up to now we have considered the simplest form of the 
two-dot spectrum, which occurs in the case when |A| is 
the largest energy scale, exceeding the tunnel coupling 
and the Zeeman splitting due to the uniform magnetic 
field by a noticeable amount. Yet it is evident that the 
qualitative features of the spectrum are sensitively de- 
pendent on the relative size of t , |A|, and Ez- As these 
parameters vary with respect to one another, the relative 
position of most energy levels can differ greatly. As a re- 
sult, there can be substantial variation in the loading and 
mixing dynamics of the two-electron states (Figs. [3]-[7]) 

In Fig. [3] we plot the two-electron spectrum of a Si 
DQD when |A| ^ i including all the Zeeman-split lev- 
els. It is interesting to follow the evolution of the two-dot 
two-electron spectrum as the ratio of the magnetic field 
to the valley splitting goes from small to large. This can 
be done by observing the way the spectrum changes from 
Figs. |3] through |6] These fi gures taken together illus- 
trate the fact that, as the magnetic field increases with 
respect to the valley coupling there is a transition be- 
tween the valley physics (which is in effect a pseudospin) 
and the Zeeman physics (which is due to the real spin). 
For large valley/Zeeman splitting the spectrum looks ef- 
fectively the same, except in one extreme case the split- 
ting between the three sets of levels is determined by | A| 
whereas in the other extreme case it is determined by 
Ez- 

On the other hand Fig. [7] contains the two-dot spec- 
trum for the opposite case when |A| ^ i. In Figs. bland 
[7] we have assumed the same value of the Zeeman split- 
ting. Finally, in general case one should be prepared for 
an intermediate scenario as in Fig. [5] in which the sets of 
levels are not clearly separated and may cross. In such a 
situation distinguishing the energy levels experimentally 
may prove challenging. 



2. Inhomogeneous magnetic fields in Si 

Quantum coherent single-qubit experiments on singlet- 
triplet qubits make use of an inhomogeneous magnetic 
field. As shown in Ref. [U an experiment on the singlet- 
triplet qubit often requires two magnetic fields: a uniform 
one to separate out the polarized triplet states, and an in- 
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0.54 0.55 0.56 0.57 0.58 0.59 0.6 

Detuning 



Figure 3: Energy level spectrum for a Si/Si02 DQD with 
a=8.2nm, b=3nm, d=2.45 and |A|=0.1meV, Ez = 0.05 meV, 
i=0.02meV. The solid lines represent singlet and unpolarized 
triplet levels S and T''\ while the dashed lines represent po- 
larized triplet levels T'^ and T~ . The top and bottom anti- 
crossings each consist of two singlets and one triplet In the 
middle anticrossing each of the three dashed lines represents 
a degenerate singlet/triplet level. 



0.4 - 




0.56 0.565 0.57 0.575 0.58 0.585 

Detuning 



Figure 4: Energy level spectrum for a Si/Si02 DQD with 
a=8.2nm, b=3nm, d=2.45 and |A|=0.01meV, Ez = 0.1 meV, 
t=0.02meV. Here the valley splitting |A| has been set as the 
smallest energy scale, the Zeeman energy as the largest scale, 
and the tunnel coupling in between. This figure illustrates 
the opposite scenario to Fig. [Sj showing the qualitatively 
different structure of the energy spectrum when |A| <t. 



homogeneous one to mix the singlet and the unpolarized 
triplet states. In the existing experiments on GaA^D the 
latter is provided by the lattice nuclear spins.'^In Si, hy- 
perfine interaction is much weaker, and number of lattice 
nuclear spins much smaller .1221 Natural Si has only 4.7% of 
the ^^Si isotope having nonzero nuclear spin (^*Si has no 
nuclear spin) and it can be further isotopically purified 
to nearly completely eliminate the nuclear spins. For the 
case when the inhomogeneous effective field arises from 




0.55 0.56 0.57 0.58 0.59 

Detuning 



Figure 5: Energy level spectrum for a Si/Si02 DQD with 
a=8.2nm, b=3nm, d=2.45 and |A|=0.05meV, Ez = 0.05 
meV, t=0.1meV. In this figure the tunnel coupling has been 
set as the largest energy scale. Although the magnitude of 
t in this graph has been exaggerated for clarity and exceeds 
what one expects to measure experimentally, this figure illus- 
trates the complications inherent in experiments seeking to 
distinguish parameters of comparable magnitude. 




L , , , , I , , , , I , , , , I , , , , tJ 

0.55 0.56 0.57 0.58 0.59 

Detuning 



Figure 6: Energy level spectrum for a Si/Si02 DQD with 
a=8.2nm, b=3nm, d=2.45 and |A|=0.01meV, Ez = 0.2 meV, 
t=0.02meV. In this figure the Zeeman energy has been set as 
the largest energy scale and the valley splitting as the lowest. 
Notice the qualitative similarity of this figure to Fig. [3] As 
the magnetic field increases, in effect the Zeeman field and 
the valley-orbit coupling trade places. 



hyperfine interaction with the nuclei in the dots, the mag- 
netic fields on the left and right dots {Bl and Bn below) 
differ by a small amount of the order of A/VN where A 
is the hyperfine constant of the material and N is the 
number of nuclei interacting appreciably with the elec- 
trons. In gated GaAs QDs this amount is on the order 
of mT. In a natural Si QD it has been estimated to be 
in the range 0.06-6 neV, corresponding to a time scale 
of O.I-lO/is and magnetic fields of 0.5-50 /iT, depending 
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Figure 7: Energy level spectrum for a Si/Si02 DQD with 
a=8.2nm, b=3nm, d=2.45 and j Aj^O.OlmeV, Ez = 0.05 
meV, i;=0.2meV. This figure illustrates the case when the tun- 
nel coupling is the largest energy scale in the problem. 



on the percentage of ^^Si in the material, from 0.01% to 
100%.!^ In systems based on purified Si, an external in- 
homogeneous magnetic field must be generated, such as 
by a nanomagnet,^" allowing better control. For exam- 
ple, one can tailor the inhomogeneous field to have only 
a f -component with a gradient along the x-direction. 

Each initialized (0,2) state has a probability of return 
to the original state after staying for a certain time in the 
(1,1) regime. In a random hyperfine field this probability 
of return is determined by the nuclear spin polarization 
difference between the two dots. Since the experiment is 
repeated many times in the time domain, the observed 
results involve averages over the inhomogeneous nuclear 
field. On the other hand, in an applied inhomogeneous 
magnetic field, the field is the same during each exper- 
imental run. Therefore, by controlling the magnitude 
and orientation of the applied inhomogeneous field one 
can control the composition of the (1,1) state after mix- 
ing, and therefore the return probability. In this case, 
if the initial state is known, we have a controlled return 
probability. The key fact is that, in the far-detuned (1,1) 
regime where (5 <C <5c, each unpolarized state (singlet or 
triplet) mixes with one other unpolarized state, while 
polarized triplets do not mix with any other state (see 
Sec. IIVB 3). Knowledge of the controlled return proba- 



bility requires exact knowledge and control of the applied 
inhomogeneous magnetic field inside each dot. 



3. Singlet-triplet mixing m an applied inhomogeneous 
magnetic field 

The inhomogeneous magnetic field is the sum of nu- 
clear and applied fields, B{r) = Bnuc{r) + Bappi{r) . Ex- 
perimentally one cannot generate an applied field with a 
profile sufficiently sharp in the f-direction to mix states 



from different valleys, so for Bappi{r) the Zeeman Hamil- 
tonian does not mix z and z states. In the case of the 
nuclear field, the contact hyperfine interaction is a sum 
of (5-functions, and can mix states from different valleys. 
Here we consider B{r) — Bappi{r). Considerations spe- 
cific to Bnucif) are discussed in the following section. 

The Zeeman Hamiltonian in Si, where g — 2 approxi- 
mately, for two electrons in a DQD is 



Hz = -f^B ki • -B(ri) + (72 ■ B{r2)] 



(18) 



where the spin operators (T12 act on the spin of particle 
1 and 2 respectively. We study the Zeeman Hamiltonian 
in the subspace of two-electron (1,1) states formed out of 
the wave functions {D±}, which is relevant to operations 
on singlet triplet qubits. Its expectation values contain 
expressions of the form (assuming weak tunnel coupling) 



{D^\Bir)\D^) = Bd 
{L^\Bir)\R^)^0. 



(19) 



We write these (real or effective) magnetic fields in the 
form Blu = {Btot ± ^B)/2, where the total mag- 
netic field is Btot ~ Bl + Br and the field difference 
is AS = B]^ — Bji. Here, since B{r) — Bappi{r), 
the magnetic field is smooth in space and has no ma- 
trix elements between states z and z. In this case, since 
{d[^1\B {ri)\Dg^^l) « 0, the matrix elements of the Hamil- 
tonian between the states {D±} are easy to find using Eq. 
pi 



{D±\Bir)\D±) = Bn 
(Z±|B(r)|i?±) «0. 



(20) 



Within the mixed- valley branch, neither an applied field 
nor the nuclear field can mix H — states with — h states 
in the mean field regime. We only consider the matrix 
elements of the Zeeman Hamiltonian within each branch 

of the two-dot spectrum, i.e. ++, and H — . 

Singlet-trip let m ixing occurs in the (1,1) configuration, 
when S ^ 111251 Without putting in valley eigenstate 
indices, we are interested in the effect of the Zeeman term 
in the space of states spanned by a (1,1) singlet "iig^ and 
a (1,1) triplet "iiif^. An applied uniform magnetic field 
separates the polarized triplet states energetically, so that 
the two-electron system remains in the branch spanned 
only by '^g^ and ^^^^f , where 



{xs\Hz\xs) = {xo\Hz\xo)=0- 
We note that 

(xsk'^^lxTo) = -z- 



(21) 



(22) 
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and the off-diagonal matrix element is 



,LR 



(23) 



The X- and y-components of AB mix the singlet with 
the polarized triplets 



LR\ 



(AS, + 2ABy) 



(vI/i«|Hz|*^f) = -^(AB.-^AB,, 



(24) 



The polarized triplets are split off by the uniform Zeeman 
field and transitions to them are suppressed. The triplets 
are mixed among themselves by the total magnetic field 



LR\ 



[B^ ± iBy 



(25) 



The uniform Zeeman field also suppresses this mixing. 
In an applied uniform magnetic field S^± mix only with 
T^i, and mix only with f^^. Each unpolarized 
(1,1) state only mixes with one other unpolarized (1,1) 
state, which is important in spin blockade. 



4-. The hyperfine interaction in Si 

The hyperfine interaction between electron and nuclear 
spins in Si is relatively well known. For donor-bound 
electrons, the full hyperfine matrix, including both the 
contact and anisotropic hyperfine couplings to the donor 
nuclei and the ^^Si nuclear spins in the envir onment, has 
been carefully measured and calculated.l^ISlHlil xhe thor- 
ough understanding of that problem allows quite precise 
knowledge of the coupled electron-nuclear-spin dynam- 
ics in Si:P, as represented by the impressive agreement 
between theory and experiment in ESR experiments on 
gj.pjso] rpj^g hyperfine interaction with conduction elec- 
trons in Si is generally known to be weak relative to con- 
fined donor electrons. The earliest experimental mea- 
surement of hyperfine coupling strength for conduction 
electrons in Si dated back to the 1950s and eOs.lS^HlSlThe 
contact hyperfine interaction, measured by average elec- 
tron probability r] = hKO)P/ (V'^) at nuclear sites, ranges 
from 178,123 to ~ 30(Pl, as compared with 77 ~ 4000 for 
GaAs."^ 

The study of hyperfine interaction in Si is currently 
drawing renewed interest because of the possibility of 
realizing spin qubits in Si. For example, in Ref,'495 
an all-electron calculation was performed to obtain the 
complete hyperfine coupling matrix (both contact and 
anisotropic parts). We note that the pseudopotential 
method is not appropriate for calculating contact hy- 
perfine interaction because it fails at the nuclear sites, 
making the use of an all-electron calculation necessary. 
The calculated rj is about 160, consistent with experimen- 
tal observations. Furthermore, the calculated anisotropic 
hyperfine coupling strength is about 3% of the contact 



part. The finite contact hyperfine coupling strength is a 
reflection that there is a finite component of atomic S- 
orbitals in the Bloch function at the bottoms of the Si 
conduction band, even though the Bloch functions are 
predominantly P-like. The small calculated anisotropic 
hyperfine coupling strength is a reflection of the highly 
symmetric nature of the Si lattice, in contrast with one- 
dimensional and two-dimensional structures such as car- 
bon nanotubes and graphene sheets. In short, hyperfine 
interaction in Si is an extensively studied subject, with its 
magnitudes in various situations already measured and 
calculated. 



An interesting contrast can be made between the hy- 
perfine interactions in silicon and low-dimensional carbon 
structures, such as graphene. We recall that the con- 
duction band minima in Si are located at ±0.85 {2'K/asi) 
along the (100), (010), and (001) axes. The six valleys are 
in the vicinity of the X points, which are high-symmetry 
points, but not at the X points themselves. Group the- 
ory tells us that the Bloch functions near the X points 
are P-like. Similarly, the Bloch functions near the con- 
duction band minima are predominantly P-like, but with 
S-like and D-like admixtures.'^ Nevertheless, group the- 
ory itself does not yield the weight of these admixtures, 
and numerical calculations are required. 

In contrast, the valleys in graphene are located at the 
high symmetry points K and K' in reciprocal space, 
where group theory tells us the Bloch function are P- 
like.l^S When this is the case the contact hyperfine in- 
teraction vanishes. The dominant hyperfine contribution 
arises from the anisotropic part of the interaction. This 
fact has interesting consequences for spin blockade, which 
were explored in Refs. 83 84, Despite the many simi- 
larities between Si and C, since the valleys in Si are not 
located at high-symmetry points, symmetry places no ad- 
ditional constraints on the form of the wave functions. In 
short, unlike in graphene, the contact hyperfine interac- 
tion in Si is finite,'^ and it actually dominates over the 
anisotropic interaction. 

For the present purpose, namely understanding the ef- 
fect of the hyperfine interaction on spin and valley dy- 
namics of singlet-triplet qubits in Si quantum dots, it is 
sufficient to think of the hyperfine coupling as the inter- 
action between an electron spin and an inhomogeneous 
Zeeman field. Both the contact and the anisotropic terms 
can be cast in this form. This field will in general differ 
between adjacent quantum dots, giving rise to singlet- 
triplet mixing when two states are close together in en- 
ergy. The hyperfine interaction can in principle mix sin- 
glets and triplets with the same valley composition as 
well as singlets and triplets with different valley compo- 
sition. These cases will be examined in what follows. We 
reiterate that any unwanted mixing present in natural Si 
can be eliminated by isotopic purification. 
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5. Hyperfine-induced valley mixing 

We write the contact hypcrfine Hamiltonian in Si as 
Hz = -Mb [o-i ■ B„y.c{ri) + (T2 ■ -B„„c(^2)], (26) 
where the hyperfine magnetic field can be written as 



Br, 



Y,^l,5{r^R,). (27) 



The contact interaction can mix states of different val- 
ley compositions. The dipolar hyperfine (anisotropic) in- 
teraction for conduction electrons in Si is two orders of 
magnitude smaller!^ so we neglect it for the following 
discussion. We again assume that a homogeneous mag- 
netic field II z that is much larger than the nuclear field 
is applied so that we can neglect mixing with polarized 
triplet states. In other words, we only consider mixing of 
S and f ° here. 

The intravalley matrix elements of the hyperfine inter- 
action have the same form as those given above for the 
applied field. We also have the intervalley term 



{D,\B(r)\D-, 



(28) 



Since we are neglecting terms of the form {L\B\R) and 
terms of linear and higher order in g and i?, {D\B\D) = 
{D\B\D) and we wiU write e.g. {b+\B\b^) = Bf_: 



Bf± ^ Bd±- {e'*"bD + e-'^"b*jj) 



1 



(29) 



In this work we take B^ and b^ as given. These quanti- 
ties can be determined by explicit evaluation of the ma- 
trix elements of Bnuc as given in Eq. (27) using the wave 



functions {-D}^ introduced in Sec. [Ill yielding 



B 



D 



AoiyQli 



i,K,Q 



fJ-B 



\F{R,rc 



K^K+Q e 



-iQRi 



>D = 



E 

i,K,Q 



Mb 



K'^K+Q ' 



^-i{Q+2kaz)-Ri 



(30) 

Beyond this level the coefficients in the expansions of 
the Bloch functions need to be determined. Very recently 
A has been worked out in Ref. 1!) which demonstrated 
that the contact hyperfine interaction is finite in a Si 
QD.I^. Explicit evaluation of b^ requires more work. 
However, inspecting Eq. (30) we expect most of the con- 



tribution to Bo to come from the term with Q — 0. 
On the other hand in bo the Q — term will be much 
smaller, due to the modulation introduced by the fast 
oscillating factor e"^'*^"-^'. The sum over i will bring this 



term very close to zero, so that we expect 6_d <C Bd. We 
stress that, since we take these quantities as given (as in- 
deed we do with A) we make no assumptions about the 
form of the Bloch functions in Si appearing in these ma- 
trix elements (except for the assumption that the contact 
hf interaction is the dominant one). 

We have the following interbranch matrix elements 



LR 



HhflT: 



LB.\ 
+ -/ 



Mb 



AB 



W++\-"hf\J- sym. 



-fJ-B 



(31) 



{Si+\Hhf\T^J}i) - ^iB 



Bi 



W 1^/1/ Ksym, 



AB 



-fJ'B 



{S"^\Hhf\T^nti) - ^Mb ■ 

In the matrix elements above the singlet and triplet la- 
bels could be reversed and the value will remain the 
same. Doing this for the first matrix element we obtain 
{T'^^\Hhf\S'^'^) = MB and likewise for the subse- 

quent ones. The hyperfine interaction does not mix the 

-f-l- and branches directly. Within the H — branch 

we have the intervalley matrix elements 

{Sl^_\Huf\Tl^) = 
(^^«|77,/|Tff ) = Bi__-Bf:^^ 

LR\ TJ \n^LR\ _ /j,LR 



(5t«ii7.^|T^f ) = {^4'\\H,M^+^) = 



1 



(S^y'LlH^flT^Jl) = ^ (AB,,++ + AS..— ) 



(32) 



syrni 2 



|F^/|Ti;« ) = - {AB,,++ + Ai3,,__) 



I qLR \ TT \rr:LR\ ^ I t}tot Tjtot \ 



The states T^^ , and S^^ , are defined in the fol- 

syinjanti sym/anti 

lowing Section. Unless specified otherwise, all remaining 
matrix elements are zero. 



C. Spin and valley blockade in multi-valley Si QDs 

Spin blockade allows charge sensors (quantum point 
contact or single-electron transistor) to differentiate the 
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two-electron singlet and triplet statesP^'i^ In this sub- 
section we explain the details of state mixing in the 

+ + / 1 ^ — branches, and the reason an inhomo- 

geneous magnetic field leads to spin blockade of only cer- 
tain states. We also consider spin blockade in a hyper- 
finc field, which has matrix elements between different 
valleys, and study their effect on singlet-triplet qubit op- 
eration. 

Spin blockade occurs when the electrons occupy the 
(1,1) state with spin/valley structure incompatible with 
the structure of available (0,2) states. In such a case upon 
changing the detuning from the (1,1) regime to the nomi- 
nally (0,2) regime tunneling of electrons into the right dot 
is suppressed, and the charge state of the DQD remains 
as (1,1). In the single-valley case this happens when we 
have a (1,1) triplet, which is lower in energy compared 
to the (0,2) triplet state after changing the interdot de- 
tuning to the detection regime, where the (0,2) singlet 
is the ground state. Therefore an initial (1,1) singlet 
state evolves to (0,2), while the (1,1) triplet state stays 
in the (1,1) configuration. This difference in the charge 
distribution can be detected by a charge sensor. In the 
multi-valley case the situation is more complicated. 

We first consider an applied inhomogeneous magnetic 

field, and study the ++/ branches, which is described 

by the Hamiltonian i?±± . For 8 the eigenstates are 
, and f^^. As the detuning is swept adiabat- 
ically from 6 ^ 6c to 6 ^ 6c, the singlets S^±. evolve 
into S^^, while the triplets f^^ remain unchanged: the 
(1,1) singlet evolves into the (0,2) singlet, while the (1,1) 
triplet remains in the (1,1) configuration, as in GaAs. 

In the H — branch, spin blockade is rather subtle 
and interesting. Consider the singlet branch [H^_ of 



Eq. (17).] The eigenstates in the (0,2) and the far- 



detuned (1,1) regimes are respectively and 



1 

75 



^sym — W + - + '^- + 1 



qLR _ J_ fqLR. _ cLR\ 
'-'anti ^ \'^-\ — — \- ) 



(33) 



In the regime of 6 ^ Sc, only S!^^ and are tunnel- 

coupled. Therefore, as the detuning is swept from 6 — 
6c <SC. to 6 — 6c ^ 0, the ground state evolves from 

into S^^, while 5^^^ remains unchanged. For the -I — 
triplet branch [H'^_ of Ec^. (17)], the eigenstates in the 



far-detuned regimes oi 6 <^ 6c and 6 ^ 6c are r|^_^ and 



rf^LR _ ^_ (rfiLR. _ rfLR\ 

anti — ^l-tj — — l) 



(34) 



As the detuning is swept from 5 <^ 6c to 6 5c: the 
anti-symmetric state T^^^ evolves into T^^ while 
remains unchanged. This is the opposite of what happens 
in the H — singlet branch. 



Operation of the singlet-triplet qubit entails initializa- 
tion of a unique (0, 2) state. When is initialized 
and the detuning is swept from 6 6c to 5 <^ 6c, this 
state evolves into S^^^. Under the action of an applied 
inhomogeneous magnetic field this state can only mix 
with the triplet f^y'^^. On the other hand, when ff:^ 
is initialized and the detuning is swept from 6 ^ 6c to 
6 6c, this state evolves into T^^^, which, in an inho- 
mogeneous magnetic field, can only mix with the triplet 
Sf^j. It is the same component of the inhomogeneous 
magnetic field that mixes S^^,^ with T^^^ and T^^^ with 

^anti- We may refer to this process as vall ey blo ckade, 
similar to that seen in carbon nanotube QDs.lSSESl As we 
sweep the detuning and allow states to mix, the destina- 
tion (1,1) states {T^^ and ^^riti) do not tunnel couple 
back to the (0,2) configuration, so that they are blocked. 
This should then allow a charge sensor measurement of 
the system states. 

Consider now an inhomogeneous nuclear field. To low- 
est order, the and branches do not mix, yet they 

can both mix with the H — branch. Some states cross, 
therefore they can be close enough in energy to be mixed. 
However the time scale for this process is on the order 
of 0.1-10 /isj^ whereas experiment uses rapid adiabatic 
passage, in which the system generally remains in the 
(1,1) regime for up to a few tens of ns, thus mixing of 
states between different branches will be avoided at the 

points where the and -I — levels cross (see Fig |3]) 

At (5 <C (5c, for mixing between different branches to be 
avoided, we need |A| 3> l-Bnud: which is satisfied by any 
A an order of magnitude larger than a few neV. 

Within the H — branch, consider initialization of a 
(0, 2) state. The singlet evolves into Sgym and the triplet 
evolves into Tanti- These two states do not mix in an 
applied inhomogeneous field, yet they can mix under the 
influence of the Overhauser field. If for example 
is initialized and evolves into Sj^^^, in the Overhauser 
field this state mixes with (matrix element Bd) or 

^aivti (niatrix element ho-) If both processes are of the 
same order, the former remains in (1,1) but the latter 
goes back to (0,2) - it is neither spin blocked nor val- 
ley blocked, and readout of the charge state of the right 
dot would not give any indication of whether mixing had 
occurred in the (1,1) configuration (yet if one starts out 
with 5*:^:^^ one can never end up with S^^ti ) We have 
seen that 6d is expected to be much smaller than B^, 
therefore the mixing of 5'f^^„ and T^^^ will take place 
on time scales much longer than 1/^s. Nevertheless, to 
determine the probabilities of return for natural Si one 
requires the intervalley matrix elements of the hy- 
perfine interaction. More work is needed to determine 
these more quantitatively. However, the ^^Si can be iso- 
topically purified, reducing the hyperfine-induced valley 
mixing as much as is needed. Under these circumstances 
we expect valley mixing to have little or no effect on spin 
blockade. 

We note that spin blockade is related to the so-called 
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supertriplet and supersinglet states in the context of 
graphene quantum dots.'^'^ In the context of Si quantum 
dots studied in this work, there are 16 states in the (1,1) 
configuration. Inspection of Fig. [3] shows that, when the 
detuning is swept from S ^ 6c to S ^ 5c, 10 of these 
states do not change to the (0,2) configuration: the nine 
triplets f^^, f^^, and f^y^ as weU as the singlet 5*^^^. 
These are the 10 supertriplet states identified in Ref. [S31 
The 6 supersinglet states, which do switch to the (0,2) 
configuration, are the singlets S^^, St-, and S^^, plus 
the three triplets T^^^- Similar conclusions regarding 
spin blockade in carbon nanotube quantum dots were 
reached in Refs. I83I84I 



D. Single-qubit gates 

For |A| ^ ksT singlet-triplet qubits in Si can be oper- 
ated in the same way as in GaAs. The physics becomes 
more interesting when the valley splitting is of the order 
of ksT (though this case makes quantum information 
processing less practical), so that any one of six initial 
states may be loaded. As the detuning is varied the sys- 
tem may go through any of the three anticrossings in 
Fig. [3) If any of the singlet or unpolarized triplet states 
is initialized, the outcome of sweeping the detuning will 
be the same as in GaAs. The polarized triplets T^^^^ 
are split off by a uniform Zeeman field that exceeds the 
inhomogeneous Zeeman field by a large amount by con- 
struction. The time scale for these states to be mixed 
with any other states by the inhomogeneous Zeeman field 
is thus extremely long, and much longer than the time 
of the experiment. As a result these states will always 
have a probability of return of 1. Overall, it is as if four 
independent qubit channels existed in a DQD, and for 
small |A| one could load any of them (with the possible 
exception of the higher-energy polarized triplet states). 
An experiment analogous to Ref. 1 is feasible, except in 
Si it will be carried out on a series of independent qubit 
channels. 

Quantum computation requires both cTj. and gates 
to achieve arbitrary one-qubit rotations. The inhomo- 
geneous magnetic field provides the ax gate, which has 
been discussed, while cr^ gates are feasible by means of 
exchange. The case | A| ^ kgT is analogous to GaAs and 
az gates can be implemented as in Ref. Jj Let us focus on 
az gates when A = 0. Once the polarized triplet states 
have been separated, four states remain (effectively) de- 
generate: S^'^ and Tfi^. If S^^ are loaded, in the 
far detuned regime they will be separated from Tj^^ by 
J. If one initializes into the H — branch one will have 
a superposition of S!^^ and T_[*_?. In the far detuned 
regime these states evolve into Ssym and Tanti- Nearest 



one can initialize into a superposition of the four states 
above. Yet each state mixes with only one other state, 
and the energy splitting between each pair of states is J. 
Thus, no matter what superposition of states is loaded, 
exchange-based Uz gates can be obtained. 



in energy to these states we have Sanu and Tsym, which 
are also degenerate, and split from Ssym and Tanti by J . 
In other words, there is only one exchange J. In general. 



It was demonstrated in Sec. Ill that exchange-based 
gates in single-spin qubits are not feasible unless |A| ^ 
ksT, whereas here they are feasible for any |A|. In 
the LDV architecture the qubits are initialized indepen- 
dently of each other. One may prepare both qubits in 
the same state (Z^ and Rj^ or and -R_), then al- 
low the exchange gate to rotate the corresponding sin- 
glets and triplets {S^± and ) into each other. If the 
two qubits are prepared in different branches we have ei- 
ther the singlets and triplets and T^:?, or S^'^ and 
T^^. We have seen above that these states are degener- 
ate: exchange splitting between is always zero. On the 
other hand, in singlet-triplet qubits, in the (1,1) config- 
uration one always accesses the states Ssym, Tsym, Santi 
and Tanti- These are superpositions of S^^ and S^^, 
and Ti^^ and Tt^. The cross terms between the H — 
and — I- states ensure that exchange between Ssym /Tsym 
and between Santi /Tanti is always J. 



E. Two-qubit operations and scalability 

One proposal for two-qubit operations has been de- 
scribed in Ref. rr31 closely related to the exchange-based 
scheme of Ref. 1 L It makes use of the fact that the 
singlet ground state is a superposition of the (1,1) and 
(0,2) configurations, while the triplet ground state within 
the operational/measurement regime is in the (1,1) con- 
figuration. The difference in charge distribution allows 
state-sensitive coupling. Coulomb interaction between 
electrons on adjacent DQDs results in an effective dipolar 
interaction that entangles qubits. In a one- valley system, 
for two qubits A and B we may construct the two-qubit 
states SaSb, TaTb, TaSb, and SaTb- If electrons 1 
and 2 are in qubit A while 3 and 4 are in qubit B, the 
matrix elements of the Coulomb interaction that lead to 
entanglement are 

Vss = {SaSb\V^^\SaSb) 

Vtt= {fAfB\Vci''\fAfB) 
Vts= {fASB\Vci''\fASB) 

VsT^ {SaTbIVJ^/ISaTb) 

Y^J" = V^ee(l,3) + Ke(l,4)+Fee(2,3) + Fee(2,4), 

(35) 

with Vec{i,j) the Coulomb interaction between electrons 
i and j. We use the generic notation Vxy for these terms, 
where X,Y = S, T. 
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In Si the DQD spectrum consists of three branches. If 
|A| < fc^T and ambiguity exists in initiahzation, then 
states from aU branches may be loaded. In Si, if both 
qubits belong to the same branch, one encounters the 
terms 

Vxy{±±.±±)^ (X±±f±±|Kn^r^s^> 

(36) 

Vxy{+-.+-)= {X+-Y+-\V,^,''\X+-Y+-)- 

Since the Coulomb interaction has no matrix el- 
ements between z and z states Vxy{++t++) — 

Vxy{ , ) = Vxy{^ — ,H — ). It does not make a 

difference which branch is loaded as long as both qubits 
are loaded in the same branch. At the same time, if the 
qubits are loaded in different branches, the interaction 
between them will be governed by matrix elements of 
the form 

(37) 

Vxy{+-.-+) = {X+-Y^+\V^^,^\X+-Y^+), 

which will all be suppressed by exponential prefactors as 
we have discussed, and will be effectively zero. The ma- 
trix elements of the Coulomb interaction between DQDs 
will be the same for qubits belonging to the same branch, 
but will be zero for qubits belonging to different branches. 
When I A| is comparable to ksT one can no longer be cer- 
tain of initializing fiducially into a particular state. The 
matrix element of the Coulomb interaction depends on 
the explicit superposition loaded. If one loads qubit A 
as aAS\'^ -I- PaS^' + 1aS\~ + S/F^^ and qubit B as 
c^bS^^ +(iBS]^ +lBS^~ +5bT^^ into the second, the 
matrix elements Vxy depend on the explicit values of the 
coefficient aA,B, Pa,b, 1a,b, 5a,b- Coulomb interaction- 
based entanglement in singlet-triplet qubits suffers from 
the same problem as single-spin qubits: it will work if 
|A| > keT, but will fail for |A| < keT. 

While single-qubit operations are formally still possi- 
ble even with the multiplicity introduced by the valleys, 
it does inhibit two-qubit operations and therefore quan- 
tum computation. The fundamental reason is that the 
Coulomb interaction between states belonging to differ- 
ent branches of the spectrum is vastly different. There- 
fore, for single-spin and singlet-triplet quantum comput- 
ers to work, it is imperative that |A| be much greater 
than kBT. We mention for completeness that the val- 
ley degree of freedom may not be a problem if two-qubit 
operations are implemented by other means, not relying 
on exchange, such as a Jaynes-Cumings type interaction 
familiar from atomic physics, yet more work is needed to 
establish a quantitative analysis of such schemes in Si. 

V. INTERFACE ROUGHNESS 

In general A is comple?f ^^ * ^°° l and is characterized by 
an amplitude and a phase. Interface roughness will in- 
troduce into the problem an additional random potential. 



Such a random potential would yield a correction to A 
and could make the amplitude and phase of A different in 
the two dots. This difference in valley composition would 
enable intervalley tunneling during interdot transitions, 
which could in principle affect two-dot dynamics, in par- 
ticular the operations of a singlet-triplet qubit. Problems 
could also arise if the exchange coupling is different in the 
three different branches -|--|-, and H — . Such a differ- 
ence could be induced by int erfac e roughness. We will 
address these issues elsewhere .1^^ Meanwhile, we would 
like to comment that interface roughness is a physical 
quantity that is difficult to determine a priori without 
destroying the device. Still, it is something that one has 
to keep in mind when considering spin physics in Si quan- 
tum dots. 



VI. OTHER TYPES OF QUBITS 

Given the additional valley degree of freedom present 
in silicon, it is natural to ask whether silicon quantum 
dots offer multiple possibilities for implementing and ini- 
tializing qubits. We discuss two other possibilities dif- 
ferent from the situation considered thus far, and show 
that, based on our current understanding and on the cur- 
rent ability to manipulate Si QDs in the laboratory, these 
alternatives cannot produce a practical qubit. 

Firstly, it is evident from our discussion that the only 

way one can be certain of initializing a particular state in 

the (0, 2) configuration without prior knowledge of valley 

splitting is to apply a strong homogeneous magnetic field 
~ I Oft 

so as to push the triplet T^_l below all the other states 
by an amount greatly exceeding ksT. The detuning can 
then be altered and one electron will tunnel to form the 
state T^l . In order to mix this state with another 
state it will be necessary to reduce the magnetic field 
greatly so that this state is close enough in energy to 
another state and the inhomogeneous magnetic field can 
mix them. In addition, the three triplet states can only 
be mixed by a uniform transverse magnetic field, but 
not by an inhomogeneous field. In the end we have a 
4-level system that has complicated and sometimes hard 
to control couplings, making it impractical for quantum 
information processing applications. 

One could also consider using the valley degree of free- 
dom to construct a qubit, in other words to employ the 
two valley-split states as a two level system and attempt 
to perform operations on it. In this case the same initial- 
ization problems described above are present. In addi- 
tion, there is practically no control over the valley split- 
ting and correspondingly no control over single qubits: 
there is no way to rotate one qubit coherently and con- 
troUably once it is initialized. The interesting paradox 
in this situation is that one could, in principle, imple- 
ment two-qubit operations by using the exchange inter- 
action provided by the valley- exchange Coulomb inte- 
gral. This valley- exchange Coulomb integral has been 
discussed above and could be used to rotate states with 
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different valley character into each other. Nevertheless, 
given the little amount of control over the valley split- 
ting, and that neither theory nor experiment have been 
able to determine the exact size of the valley-exchange 
Coulomb integral jy, this would be a highly impractical 
method at present. 



VII. IDENTIFYING VALLEY-SPLIT STATES 

To implement a qubit in a reliable manner we need to 
determine whether the condition |A| ^ ksT is satisfied. 
We recall that A has been calculated recentljP^ and it 
is believed that its magnitude can be as high as «0.25 
meV, and can be further increased by applying an electric 
field. In Ref. [52| we described a quantum coherent ex- 
periment to measure the valley splitting. An alternative 
approach to measure |A| is to use a two-electron single 
quantum dot and probe its ground state while sweeping 
the magnetic field. The basic idea here again relies on the 
detection of the magnetic field at which the ground state 
undergoes a singlet-triplet transition. This single-dot ex- 
periment requires charge sensing. One should adjust the 
gate potential in such a way that the dot contains two 
electrons, so that we can refer once more to the spectrum 
of the doubly-occupied single dot as in Fig. [2j One can 
then perform tunnelling spectroscopy to determine the 
spectrum of the QD, and identify the ground state. At 
zero applied field the ground state should be a S sin- 
glet. As the applied field is increased, the ground state 
remains to be the singlet and barely changes in energy 
with the field, until it switches to beyond which 
the ground state energy energy would decrease linearly 
with the increasing field. This singlet-triplet transition 
in the ground state occurs when Ez = 2|A|. Recording 
the B-ficld dependence of the ground state energy and 
identifying B for the singlet-triplet transition gives a re- 
liable estimate of |A| in this QD. This description also 
shows that the experiment can be reformulated alterna- 
tively as a transport (resonant tunnelling) experiment. 
Recent experimental work on Si QDs^ suggests that 
this experiment is feasible at present. We are indebted 
to discussions with Malcolm Carroll for this idea. 



VIII. SUMMARY AND CONCLUSIONS 

We have analyzed the way the valley degree of freedom 
affects spin qubits in Si quantum dots. The multiplic- 
ity of the ground state alters the energy level spectrum 
that is accessed in the implementation of single-spin and 
singlet-triplet qubits. Moreover, unlike in single-valley 
systems, the energy scale is no longer set by the con- 
finement energy, but by the valley splitting, which is ex- 
pected to be approximately one order of magnitude lower. 

Initialization is a crucial step for any qubit. Fo Si QDs, 
when |A| ^ ksT, the system is effectively a one-valley 



system, and both single-spin and singlet-triplet qubits 
are analogous to those in GaAs. 

When the valley splitting |A| is on the order of fc^T 
or smaller, both single-spin and two-spin qubit architec- 
tures become essentially impossible. In particular, for 
a single-spin qubit, different single-electron valley eigen- 
states cannot be energetically distinguished, so that reli- 
able orbital initialization cannot be realized even if spin 
states can be properly prepared. While single-spin rota- 
tions can still be performed, the existence of non-identical 
orbital states will make the exchange coupling between 
neighboring quantum dots unpredictable and dramati- 
cally increase the difficulty of two-qubit gates. Therefore 
spin quantum computation becomes practically impossi- 
ble in this limit. 

In singlet-triplet qubits, if |A| < ksT and provided 
a uniform Zeeman field separates the polarized triplet 
states, any one of four two-electron states may be loaded 
on a single dot: two singlet states in which the electrons 
are in the same valley eigenstate, as well as one singlet 
and one unpolarized triplet state in which the electrons 
are in different valley eigenstates. The DQD spectrum is 
characterized by three distinct anticrossings, which rep- 
resent three different branches of the spectrum, each with 
a different valley composition. Depending on the relative 
values of |A|, the Zeeman energy Ez, and tunnel cou- 
pling t, these branches may interpenetrate, resulting in a 
complex experimental energy level spectrum. An applied 
magnetic field does not mix valley eigenstates with differ- 
ent valley composition. In a nuclear magnetic field the 
dominant term is due to the contact hyperfine interac- 
tion. An intervalley term exists that in principle lifts the 
spin blockade in the middle (valley-mixing) branch of the 
spectrum. We expect this term to be much smaller than 
the intravalley matrix element of the hyperfine interac- 
tion, which is the relevant matrix element for gates. 
Therefore, although singlets and triplets from different 
branches of the spectrum may mix in an Overhauser field, 
this mixing occurs on time scales much longer than those 
relevant to quantum coherent experiments and ax gates 
are feasible. Furthermore, since exchange is the same in 
all three branches of the DQD spectrum, CTz gates can also 
be implemented even if |A| < fc^T. For |A| comparable 
to ksT the system could be viewed as a set of four inde- 
pendent qubit channels. At the same time, the strategy 
for scaling up singlet-triplet qubits proposed in Ref. [T3| 
is not practicable when |A| < fcsT, since the Coulomb 
interaction is again unpredictable between neighboring 
qubits. Yet spin blockade in Si in an applied inhomoge- 
neous magnetic field occurs in the same way as in GaAs. 

For both single-spin and singlet-triplet qubits, one- 
qubit operations are feasible in the presence of a val- 
ley degree of freedom. Two-qubit operations based on 
the Coulomb interaction are only possible if both qubits 
have the same valley composition, and depend on the 
ability to initialize all qubits unambiguously into the 
same state. Architectures in which entanglement relies 
on on the Coulomb interaction are not scalable unless 



16 



|A| ^ ksT. For these schemes the only feasible option is 
to ensure |A| 3> ksT and implement spin qubits in the 
same way as in a single-valley system such as GaAs. 

We have shown also that other types of qubits involv- 
ing the valley degree of freedom are not feasible based on 
our current knowledge and ability to control Si QDs in 
the laboratory. 

Valley-split states and the magnitude of the valley 
splitting may be identified by sweeping a magnetic field. 
This can be done using charge sensing or transport of 
two electrons on a single quantum dot, or by performing 
a quantum coherent experiment on a singlet-triplet qubit 
in a double quantum dot.l^ These measurements can be 
checked against one another. 

The valley degree of freedom of Si poses different qual- 
itative problems for Si QD spin qubits compared with 
the corresponding Si:P bulk spin qubits, where valley de- 
generacy leads to the exchange oscillation phenomenon, 
making the fabric ation a nd control of two-qubit exchange 
gates problematic! ^^ l ^°^ In that case, it is the lifting of the 
valley degeneracy (rather than the valley degeneracy it- 
self) by the singular Coulomb potential of the donors that 
is at fault. It appears that the valley degree of freedom is 
a serious enemy of solid state spin quantum computation 
independent of whether one considers surface or bulk spin 
qubits. In fact, other solid state systems where valley de- 
generacy exists, such as graphene or carbon nanotubes, 
may also have difficulties with respect to spin qubit ar- 
chitecture because of the competition between spin and 
valley degrees of freedom. 



For a few qubits the existence of valleys closely spaced 
in energy is not necessarily a problem, since experiment 
can use trial and error and explicitly characterize each 
sample, ensuring that the valley splitting is sufficiently 
large. Nevertheless the quantum computer architectures 
discussed in this work are not scalable unless one de- 
termines a reliable way to control the valley splitting 
systematically. The ideal solution would be to develop 
a quantitatively precise control over valley splitting so 
that it can routinely be engineered to be large. Unfortu- 
nately, designing samples with large valley splitting thus 
remains an elusive task. 
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Appendix A: Matrix elements 



The interdot tunneling matrix element t has a single-particle part to and a Coulomb-interaction (enhancement) part 
s, and we write it as t = io + s. Starting with the wave functions \D^) for a single electron in a Si DQD, the matrix 
elements of the Hamiltonian within the Hund-MuUiken approximation involve the following integrals (the superscript 
(i) labels the i-th electron) 



to = {L^\Ho\R^} 

k = (if ) = {LfR^^l\V..^Lf^R^^l) 

J (i«i?f iKelif i?^^) = {Lf^R^^l\V.e\L'^lRf^) (Al) 
. = (Lf^LfW^^lLf^Rf) = (LfL^^lW^^lLf^R^^l) 



We shall refer to these as the bare matrix elements. On an asymmetric double dot we have to allow for the energies 
on the left and right dots, El and er, to be different, as explained in the text. When one switches from {\L), \R)} to 
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the orthogonal basis spanned by {|Z/), the matrix elements of the Hamiltonian involve the following integrals 

_ _ £l + g^SR ~ 2gto 
l-2lg + g^ 

. ^ gfl + g^£L - 2gto 
1-2/5 + 52 

r _ tQ{l + g^)-g{£L+£R) 
l-2lg + g^ 

- _ fc(l + ff^) - Agsjl + g') + Vi + 2g^u 

(l-2/.g + .g2)2 

^ ^ i(l + g^) - 4g5(l + g^) + 2g2(fc + j + u) 
3 ii_2lg + g^Y 

. ^ a(l + 6g2 + g^) - g(l + g2)(2j + fc + ^) 
^ (l_2/g + g2)2 

. u{l + - 4y.s(l + //) + 4y-',y + 2r/(A- + 2j) 
(l-2/g + g2)2 



Appendix B: Coulomb integrals 

Most generally, we wish to evaluate the integral 

Define 

F(fe) = / S r f(r). 

The Fourier transform of /(r) = 1/r is -F(fe) = A-K/k"^. We write the term l/|ri — as a Fourier expansion 

A-TT r p-ik2{ri-r2) 1 r p-ik2-(ri-r2) 

'd^k,^—^ = ^Jd^k,'-^. ■ (B3) 



In -rsl (2^)3 J ' kl 27r2 7 ' kl 

Next we Fourier transform all the other terms (Gaussians) entering /(P, Q). First we take the 2;-dependent terms 

(B4) 

2V7r J 

We also Fourier transform the x and y-dependent terms 



„2. 2 



(B5) 



The integral I{P, Q) becomes 



d^rid^r2d^kid^k2d^k3 e-'('=i+'=^)"-ie-'('=3-'=^)"'= p_:/c. o_ (°^^L+'-^'°L) .(^Mx+i^ 

(B6) 



1 M!fci^ife,.(P-Q) -<^^M4i^ 



27r2 y /fc2 
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The integral depends only on the difference R = Q — P. If i? is in the xy-pl&ne we obtain 



C""-"")''! _ „ „ fbk 



dk±e-" — ^^Joik±R)Erfci^^j. (B7) 

In 2D, i.e. as 6 0, wc have /2d(0) = ^-v/f^- In Hund-MuUikcn problems wc typically require evaluation of four bare 
Coulomb integrals (direct and exchange, on-site repulsion, plus an additional term), namely 

2 2 

k=— I{-Xox,Xox) = — I{2Xox) 



2 2 



(B8) 



w = — I{XoX, Xqx) = — /(O) 

2 2 

s = ^ Z I {-Xqx, 0) = ^ / /(Xox). 

Appendix C: Matrix elements in a magnetic field 

In a magnetic field the wave function of each dot acquires an additional phase, and we write them as L±m, R±m 

(CI) 

The magnetic length is = ^/K/eB, the Larmor frequency = eB /2m, the effective frequency lo = -^/wq + and 
the parameter hm = ^/l + uj\/ijJq, while a, d and I all change to am = o-l\fWrn d,n = d\/h/a, and Im = (26m-i/i'm)_ 
The Coulomb integral I changes to Im, in which all the parameters are replaced by their corresponding expressions 
in a magnetic field, so a ^ and so forth. The Coulomb matrix elements change as follows 

km = — Imi^Xox) 



,2 T ( ^^0<^m -X 



_ '2 

e 



(C2) 

e 

Sm = —Im Im{XoX —3 y)). 

6 R 



The bare tunneling parameter to changes to e tm, with found from to by replacing a with and so forth. 

Appendix D: Valley-exchange Coulomb integral 

The valley-exchange Coulomb integral for a dot located at Rd =0 

= J d'n J d\, D^'^Df^Vee Df^Df). (Dl) 
The aim is to reduce this integral to the form of expression xI{P, Q) as calculated above. Firstly, 



^ vee-f - _.-^„4i.2 ZiZ ITT- - c - ^ul{ri)u^{ri)ul{r2)u^{r2). (D2) 
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We expand the lattice-periodic wave functions as 



















ulir2) 






K2 
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where /«i = Ki — K3 and K2 ~ K2 — K^. Abbreviating ^ = X^jCi ff2 k,i 1^1 + "^kz — Hi and k,2 ~ K2 — 2kz 
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(D3) 



The integrand is almost in the form we want it, except we need to complete the square in the Gaussians. To shorten 
the algebra a little we note that the only sizable terms in the sum over reciprocal lattice wave vectors will have 
Kix = Kiy = ^2x = = 0. All the rest of the terms are exponentially suppressed. Consequently we only need 
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Therefore the valley exchange Coulomb integral can be expressed as 
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Using the result for the integral / found above and Akz — Riz — R2z 
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We redefine kz ^ kz + in the above integral without affecting the limits of integration 
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where afc_L = q_L, bkz = Qz and bAnz = Aqz, so that 
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(D8) 



^2/ I ^ \2 

The exponent in the overall prefactor reduces to i.'^i^+"^2W ^ ^j^g prefactor has no dependence on feg. The 

coefficients for Si were calculated in Ref. £8_ using a pseudopotential method. For the final step in the evaluation 
of the integral we examine which reciprocal lattice vectors make a sizable contribution. We need + K2z = 0. The 
first possibility is = k.2z = 0, with Aqz = 46fco, which yields, using the coefficients found in Ref. 1981 



E ~ E ^KiCKyK2'^K2 
Ki,K2 



0.01. 



(D9) 
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The second possibility is ki^ = — K2z, in which case Aq^ = 2bKiz +4&fco. These terms cancel exactly because of lattice 
symmetry. For example, consider the case kiz = 2 {271 /asi) (that is, the Umklapp terms), for which 

Umklapp Ki,K2 

Using the fact that cZx = c|f , the sum of the largest coefficients entering J2umklapp zero. 
The valley-exchange Coulomb integral is approximately 



3v^ — 
ea 



(E) -/ dq^<l±^~^ dq. ^ (Dll) 



For a = 8.2 nm, b = 3nm and e « 7.9, we get Ag^ « 117 and (e^/ea) « 22.2meV. The integral in brackets is evaluated 
using Mathematica and yields « 8.5 x 10~^, so the final res ult is jy » 0.02/xeV. 
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